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Abstract 

In this short article we discuss the properties of electromagnetic fields in heavy- 
ion collisions and consequences for observables. We address quantitatively the 
issue of the magnetic field lifetime in a collision including the electric and chiral 
magnetic conductivities. We show that for reasonable parameters, the magnetic 
field created by spectators in a collision is not modified by the presence of 
matter. Based on this, we draw definite conclusions on observed effects in 
heavy-ion collisions. 



1. Introduction 

The experiments with heavy-ion collisions of ultra-relativistic energies probe 
not only matter in extreme temperatures and densities, but also under action of 
extremely strong electro-magnetic fields with magnitude of the hadronic scale, 

eBr^mi [ma. 

The magnetic field is a key ingredient for many observables related to local 
parity and charge parity violation [1^ . The lifetime of the magnetic field, which is 
needed to describe the observed data of the elliptic flow dependence for positive 
and negative charged particles on the asymmetry, within the framework of the 
Chiral Magnetic Wave must be as large as 4 fm/c [3]. However, the photon 
azimuthal anisotropy measured at the top RHIC energy can be described by the 
magnetic field at a time scale of a few 0.1 fm/c This apparent discrepancy 
demands theoretical studies of the time dependence of the magnetic field. In 
this short article we will consider effects of finite conductivity on the lifetime. 
We also discuss the dependence of the magnetic field on the collision energy 
and draw some conclusion on photon azimuthal anisotropy for RHIC and LHC 
energies. 
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2. The Electro-Magnetic field in heavy-ion collision 

The Maxwell equation describing the time evolution of the electromagnetic 
field in a collision reads 

f = -VxE, (1) 

^=VXE-J; (2) 
where the electromagnetic current can be decomposed in two pieces: 

j=jcxt+jint, (3) 

the internal current, jint, and the external one of the colliding nuclei, jext- 
The latter we will treat in the eikonal approximation neglecting eff'ects of the 

proton deceleration and/or stopping. For later convenience let us separate the 
electromagnetic field also in the "external" and "internal" , e.g. 

E = i?ext + -Eint, (4) 

where the external electromagnetic field satisfies the equations 



dt 

dt 



while for the internal we get 



dt 

dt 



-V X i;ext, (5) 

VxBext-Jext; (6) 

■ -V X 4nt, (7) 

V X Ant -Jint; (8) 



This representation is especially convenient for finding a numerical solution due 
to the following two reasons. First, there is no need in solving the first couple of 
equations for the external components with the singular source terms, since the 
solution can be obtained by boosting the electric field of both nuclcnis. Second, 
the singularities of the sources are spread by the fields, which leads to a better 
convergence of a numerical scheme. 

In the matter rest frame, the internal current may include the contribution 
from the Ohmic conductivity 

Johm = cr£ (9) 
and the one induced by the QED anomaly 

ianom = CT^-B, (10) 



2 



where tr^ oc /is and is the axial chemical potential, that can be created by 
a sphaleron transition in a collision. The value of is quite uncertain, so is 
a^. The electric conductivity crohm was calculated using a first principal lattice 
QCD approach in the quenched approximation. It has been found |S] that above 
the transition temperature 

<^o?m° = (5-8±2-9)|^MeV. (11) 

c 

The quenched approximation used in Ref. _5j does not allow to draw a final 
concussion on the magnitude of the conductivity. We, however, do not expect 
that taking into account the quark loops will modify the results larger then by 
a factor of 2. 

We also assume that charged matter, if created at early times, does not 
develop large collective velocity v, which allows us to neglect the contribution 
V X B when transforming to the laboratory rest frame. Both conductivity iJohm 
and may result in a substantial increase of the lifetime of the magnetic field. 
In what follows we consider both cases separately. 

Magneto-hydro scenario, — and crohm ^ where is the charac- 
teristic timescales. For this case, the Maxwell equations can be written as 

V Bint - + ^'^^"'~df~ + ^ ' 

For late times the external contribution, the last term, is not important and can 
be skipped. The first term on the right-hand side of Eq. ( 12 ) can be neglected 
for (Tohm 3> 1/tc yielding a diffusion equation for the components of the inter- 
nal magnetic field. The diffusion equation leads to a significant increase of the 
lifetime of the magnetic field, as was pointed out in Ref. [7] recently. However, 
this result was obtained under the assumption of crohm ^ l/^c- The character- 
istic time scale is defined by the external magnetic field and proportional to the 
sickness of the nucleus in the beam direction, i.e. tc ~ 2R/^. For the top RHIC 
energy, tc ^ 1/5 fm = 10"^ MeV"^ The condition aohm > l/^c ~ 10^ MeV, 
therefore can be rewritten as a condition on the temperature of the created 
plasma, T/Tc ^ 10^, which is unrealistic for collisions of such energy. Thus we 
conclude that for realistic values of the parameters the condition aohm 3> 1 /tc 
is never satisfied. In the next section, we will demonstrate our numerical results 
for realistic values of crohm) illustrating the same fact. 

Formation of knots, ay, ^ 1/tc and (Johm = 0. In this case, neglecting second 
derivatives with respect to time we get 

V^^int + Cr^V X Bint + fT^V X Bext = (13) 

As in previous case, for late times, the external contribution can be neglected 
resulting in 

Sint = V X Bint. (14) 



Similar result was obtained in Ref. [6] earlier. 
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Figure 1: Magnetic field for static medium witli Ofimic conductivity, cohn 



The total magnetic helicity 



H 



d^xA B 



(15) 



is conserved for the closed systems. For two fluxes of the magnetic field (pi and 
<j)2j the helicity can be related to the linking number H — 2n(j)i(f)2- Substituting 



helicity 



Eq. (14) to Eq. (151 and performing trivial transformations we obtain for the 
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H^ = / B^d^x = Eb, (16) 

f^x Jv o-x 

where Eg is the total magnetic energy. This shows that conservation of the 
helicity leads to the conservation of the total magnetic energy for the processes 
with the timescales, tc <^ l/^x- The volume V in Eq. (16 1 is defined by the 
region of space, where cr^ ^ 0. Owing to the expansion of the medium this 
volume grows in time roughly as for late times, as t for early times. Therefore 
we expect the magnetic field to decay according to the power law B ~ 
or i? . This is somewhat slower then the decay of the field induced 



by the spectators B. 



spcct 



r 



This conclusion, however, does not take into 



account the formation of non-trivial topological objects, knots of the magnetic 
field with non-trivial linking number. As was shown in Ref. j8], the higher 
the linking number corresponds to longer lifetime of the magnetic field up to 
Bn ~ t~^^^. Returning to the constraint cr^ ^ l/^ci we can roughly estimate 
if this is satisfied in heavy ion collisions. The chiral conductivity is defined by 
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chiral chemical potential following [3] 



-x=(^^^cE'?/jM- (17) 

To make an optimistic estimate we will consider quite high values of /i5 ~ 1 
GeV. Nevertheless even in this case the numerical value for is only 15 McV, 
while the inverse characteristic time for a collision of the top RHIC energy is 
10'^ MeV. This again shows that the effects of finite cr-^ will not be important 
for the top RHIC and LHC energies. 



3. Numerical results 

Numerical results to be presented in this section are obtained with the Yee 
algorithm ,10j . which is numerically stable for the conductivity ranging from 
to CO. The calculations proceed as follows. First, we initialize the distribution of 
protons both in projectile and target according to the Woods-Saxon distribution 
with the standard parameters TT]. We also checked that up to statistical fluc- 
tuations the results obtained using the Wood-Saxon distribution coincide with 
those in approximation of collision of homogeneously charged spheres, which 
will be used for the sake of illustration in Fig. 1. Next we assume that both 
target and projectile move with the opposite velocities of the same magnitude 
v"^ = 1 — (2mp/-y/s)2. The conducting medium in the collision is not formed 
immediately, because the quarks need time to be created from the glasma field. 
Nonetheless, to make our estimates of the conductivity effects as optimistic as 
possible we will consider that the conducting medium is formed immediately 
after the collision and does not alter during the evolution. We also neglect 
possible non-equilibrium effects and finite time response of the medium to the 
electric field, which as was shown in Ref. '12 may play an important role. 

The decay of the conductivity owing to expansion of the medium can only 
decrease the lifetime of the magnetic field and thus will not be considered here. 
Our simulations are done for Au-Au collisions at energy ^/s = 200 GeV and fixed 
impact parameter b — Q fm. In Fig. 1 we show time evolution of the magnetic 
field in the origin x — Q a& a, function of the electric conductivity crohm- The 
results show that the lifetime of the strong magnetic field {eB > m^) is not 
affected by the conductivity, if one uses realistic values obtained in Ref. [5j . 



4. Energy dependence 

In the previous section, we established that for realistic values of the conduc- 
tivities the electromagnetic fields in heavy-ion collisions are almost unmodified 
by the presence of the medium. Thus one can safely use the magnetic field gen- 
erated by the original protons only. This magnetic field can be approximated 
as follows 
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Figure 2: Coupling of gluon field to photons. The cross labels external coherent magnetic 
field line. 

where Z is the number of protons, i? is the radius of the nuclei, 7 is the Lorentz 
factor and, finally, c is some non-important numerical coefficient. We are inter- 
ested on the effect of the magnetic field on the matter. Thus we need to compute 
the magnetic field at the time tm, characterising matter formation time. On the 
basis of a very general argument, one would expect that t„i = aQj^. The phe- 
nomenological constraints from photon azimuthal anisotropy at the top RHIC 
energy demand ~ 2i?/7RHic, i-e. a — 2i?Q^^^'-'/7RHic- Using this relation, 
we can estimate the magnitude of the magnetic field at the LHC energies at the 
time t — tjn- For the LHC, since t„i <C 2i?/7 we have 



Using Qs oc (v^)°'^, we get (ei3LHc/ei?RHic)^ ~ 0.32. Therefore the magnetic 
field is not suppressed strongly at matter formation time at LHC. However, 
the anisotropy of photon production from magnetic field is suppressed. This 
follows from the following dimensional argument. First of all, to couple gluons to 
photons we are considering the box diagram showed in Fig. [2] (see also Ref. [i]). 
Integrating quark fields, we obtain an effective coupling g'^e^hG^F^ , where = 
F^^F^^ and similar for . At early times, according to saturation picture the 
contribution, g^G"^ is of the order 1 in powers of .g^. One of the factors of F^^, is 
associated with external coherent magnetic field, thus it compensates one power 
of e. Therefore, the diagram is of the order g°e, i.e. of the same as Coulomb or 
Compton contributions to photon production. The prefactor h is defined by the 
momentum running in the loop. In the first approximation it is proportional to 
Q^"* and canceled out by the gluon G^ . Finally, parametrically we have that 
the amplitude of photon production is eBk± (we consider only central rapidity 
region, kz =0). The rate of photon production is thus [eB^k^. 

Thus the photon V2 generated at early stage is thus, B'^k\/{Q^g). Therefore, 
the ratio 




(19) 




(20) 
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calculated at given k± is small owing to the prefactor ( qlhc j ■ 

5. Conclusions 

In this article, we calculated the time dependence of the magnetic field on 
time allowing for conductivity effects in the plasma. In contrast to the results 
obtained in Ref . [7] , we showed that the effects of conductivity do not play an 
important role for realistic values. We also have argued that while photon flow 
may receive reasonable size effects from production in an external magnetic field 
at RHIC energies, the induced flow is very small at LHC energies. This can also 
be tested experimentally with the method proposed in Ref. [13]. 
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